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Abstract 



We study classical and quantum dynamics of two spinless particles confined in a quantum wire 
with repulsive or attractive Coulomb interaction. The interaction induces irregular dynamics in 
classical mechanics, which reflects on the quantum properties of the system in the energy level 
statistics (the signatures of quantum chaos). We investigate especially closer correspondence be- 
tween the classical and quantum chaos. The present classical dynamics has some scaling property, 
which the quantum counterpart does not have. However, we demonstrate that the energy level 
statistics implies the existence of the corresponding scaling property even in the quantum system. 
Instead of ordinary maximum Lyapunov exponent (MLE), we introduce a novel kind of MLE, 
which is shown to be suitable measure of chaotic irregularity for the present classical system. We 
show that tendency of the energy dependence of the Brody parameter, which characterizes the 
energy level statistics in the quantum system, is consistent with that of the novel kind of MLE. 

PACS numbers: 05.45.Pq, 45.50 Jf, 73.21.Hb, 73.23.-b 



I. 



INTRODUCTION 



* syunpei710@cmpt.phys.tohoku.ac.jp 



The recent development in high technol- 
ogy has fabricated nano-scale quantum dots 
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containing a finite number of interacting elec- 
trons and optically trapped atoms where a 
finite number of interacting macroscopic par- 
ticles are trapped in a small area. Quantum 
mechanics of these systems constitutes a top- 
ical subject. In such systems, the underlying 
classical motion is expected to play an im- 
portant role. The nature of the classical mo- 
tion, i.e. regular, mixed, or chaotic charac- 
ter, reflects on some of the quantum proper- 
ties of the systems, particularly in the energy 
level statistics. In this context, a number of 
studies on the quantum chaos of systems con- 
taining a few electrons in quantum dots have 
been reported pj{9j. However, closer corre- 
spondence between the classical chaos and 
the quantum chaos in those systems has not 
been investigated well. 

The simplest system among them would 
be the one-dimensional system J5j|6]. In this 
paper we are concerned with behavior of two 
particles interacting with each other via the 
repulsive or attractive Coulomb potential in 
a one-dimensional system and study the cor- 
respondence between the classical and quan- 
tum chaos in detail. According to Fendrik et. 
al. [6], we introduce an effective Hamiltonian 
for a quantum wire, which reduces the orig- 
inal 3D system to the quasi-one-dimensional 
system. While its classical dynamics has 
some scaling property, the quantum coun- 
terpart has no such scaling property. We, 

2 



however, show that the energy level statis- 
tics implies the existence of the correspond- 
ing scaling property even in the quantum sys- 
tem. This is demonstrated by calculations 
of the Brody parameter for distributions of 
the nearest neighbor level spacing (NNLS). 
This subject, scaling in quantum chaos, has 
been examined for some other systems, cou- 
pled harmonic or quartic oscillators [T0l - [T2] 
and the hydrogen atom in a magnetic field 

In order to clarify closer correspondence 
between the classical and quantum chaos, 
we introduce a novel kind of maximum Lya- 
punov exponent (MLE) instead of ordinary 
MLE. The ordinary MLE is a measure of the 
rate per unit of time for separation between 
two adjacent orbits while the new MLE is 
the one per unit of distance for separation be- 
tween them. The new MLE is a suitable mea- 
sure to compare chaotic irregularity among 
classical orbits with different energies. We 
show that tendency of the energy dependence 
of the Brody parameter is consistent with 
that of the new MLE. We further show that 
the area of a chaotic region in Poincare maps 
are not a suitable measure of chaotic irreg- 
ularity for the present system, while several 
authors showed that it is a suitable measure 
of the irregularity in other systems [T3l - [T5] . 

This paper is organized as follows: In 
SecjTTJ we construct a quasi-one-dimensional 



model of two electrons confined in a quan- system is written as 

turn wire. We introduce a new kind of ^-^ r h 2 f d 2 d 2 d 2 

MLE. In SecjTTTl we explore the distribution ^ 2 L 2m \dx 2 dy 2 dz 2 



of NNLS in wide range of energy and interac- _^I mCc , 2 ( x 2 _|_ y 2 ^ 

tion strength. Then we examine the chaotic g 2 



±- 



irregularity of the corresponding classical sys- ~ x i) 2 + (Vi ~ V2) 2 + (>i - z 2 ) 2 ' 

tern with the use of the MLE and Poincare (1) 



maps. We clarify correspondence between 
the energy dependence of the distribution of 
NNLS and the chaotic irregularity in the clas- 
sical counterpart. Summary and conclusion 



are given in Sec, IV 



We assume that the particles occupy the 
lowest-energy state associated with the trans- 
verse motion, which is energetically well sep- 
arated from the excited states. Then the two- 
particle wave function can be approximated 
as 



II. MODEL AND METHOD 



A. Quantum dynamics (2) 

We consider two spinless particles (two where (j) (x) is the lowest energy eigenstate 

electrons or an electron-hole pair with the of a harmonic oscillator. The wave function 

same mass) confined in a quantum wire. We $(^1, z 2 ) satisfies the equation 
assume a narrow parabolic confinement in the 

H 1D <f>(z u z 2 )=E<f>(z u z 2 ) : (3) 

transversal directions (x and ^/-directions), 

which are much narrower than a confinement where the effective Hamiltonian H 1D is de- 

in the longitudinal direction (z-direction). fined by 

We consider a hard wall potential in z- TT h 2 ( d 2 d 2 \ Tr 

HlD = la? + a? + Vid(1zi ~ Z2iy 

direction. The particles are interacting with V^i u ^2/ 

(4) 

each other via the repulsive or attractive 

Coulomb potential. The Hamiltonian of the Vid(z) is the effective potential given by 
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V 1D (z) = ±e 2 



= ± 6 - 
a 



\Mxi)\ 2 \Myi)\ 2 \M^)\ 2 \My2)? 

yj (xi ~ x 2 ) 2 + (yi - V2) 2 + z 2 
00 sexp[-s 2 /2] , 

— CIS 

^s 2 + ( z / a ) 2 



dxxdyidx 2 dy 2 



(5) 



and a = \jYijmuj. In this way, our sys- 
tem exhibits a quasi-one-dimensional prop- 
erty. Now we introduce a model potential [6] 
defined by 



V m (z) = ± 



(6) 



In Fig{TJ the solid and broken curves indicate 
the numerically calculated potential V\d and 
the analytical potential V m) respectively. We 
see that V\d can be well approximated by V m . 
We adopt V m instead of Vm as the interaction 
potential between particles since the analyt- 
ical potential can be dealt with more easily. 
Thus our effective Hamiltonian is written as 



eff 



2m\dzf dz\ 

±, e2 . 

^a 2 + <yZl - z 2 ) 2 



(7) 




We scale lengths, angular momentums and 
masses by L, h and 2m, respectively, where 
L is a width of the system in the longitudinal 
direction. Then the effective Hamiltonian is 
reduced to 

A 



H. 



d 2 d 2 

+ 7T^7 I + 



e// 1 dz 2 ' dz 2 J ' - Z2 ) 2 + S 2 

(8) 



FIG. 1. The numerically calculated potential 
Vid (solid curve) and the analytical potential 
V m (broken curve). The lengths and potentials 
are scaled by a and e 2 /a, respectively. 

where A is the effective interaction strength 
parameter given by A = ±2e 2 mL/h 2 and 
S = a/L. Note that the parameter A can be 
expressed as A = ±2L/a#, where a# is the 
Bohr radius defined by a# = h 2 /e 2 m. The 
particles are confined by hard walls within 
< z x < 1 and < z 2 < 1. These hard 
walls describe the boundary of the quasi-one- 
dimensional wire. We examine the feature 
of the system as varying A while keeping 5, 
which implies that we vary the system size 
keeping the ratio between the longitudinal 



and transversal lengths. the cases of fermions. 

In order to look for signatures of quan- 
tum chaos in the present system, we exam- 
It should be noted that the present ID ine distributions of the nearest neighbor level 
two-particle system is equivalent to a 2D spacing (NNLS). The eigenenergies are ob- 
system of one particle having a coordinate tained by diagonalizing the Hamiltonian ma- 
(zi,z 2 ) within a hard- walled square billiard, trices numerically, whose elements are evalu- 
The Hamiltonian of the latter system is also ated by using energy eigenstates without the 
given by Eq.([8]), in which the first and sec- Coulomb interaction (Slater determinants) as 
ond terms represent the kinetic energy of the a basis set: 

particle and the third term represents a ex- / 

0m,n(^ij Z2) = v2( sin(m7TZi) sin(n7T2:2) 
ternal potential. We can chose energy eigen- ^ 

functions of the 2D system as being sym- sin(ri7rzi) sin(m7rz2)^ , (9) 

metric or antisymmetric against exchange be- where m and n are integer larger than zero, 
tween z x and z 2 , i.e., ^(z u z 2 ) = ^(z 2 ,z 1 ) The components of the Hamiltonian matrices 
or $(^1,^2) = ~^(z 2 ,Z\) ) since the Hamil- are represented with respect to m n in Eq.g 
tonian does not change under this exchange. as 

The symmetric and antisymmetric cases cor- 2/2 2\ e 

< (p m ,n\H\(p m >y >= 7T (m +n)d m ^d n y 

respond to the boson and fermion cases, re- 

+4A{/(m, n\m f ) n') — /(m, n\n' ) m')}, (10) 

spectively, in the ID two-particle system. In 

the present paper we are concerned only with where /(m, n\m' , n') is defined by 
1 



I(m,n\m f \n f ) = / / dz\dz 2 
Jo Jo 



sin(m7TZi) sin(n7rz 2 ) s'm(m f 7rzi) sin(n / 7rz 2 ) 



Jo y/{ Zl ~ Z 2 y + (52 

(11) 



We further take into account the parity of (odd, odd) and those having odd parity with 
the system. The present system is invari- (m,n) =(even, odd) or (odd, even). We con- 
ant under the inversion associated with the centrate ourselves on the eigenstates of the 
center (zi,z 2 ) = (1/2,1/2). Therefore, the even parity in this paper, when we examine 
eigenstates are classified into ones having the NNLS. 
even parity with (m, n) =(even, even) or 
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NNLS is fitted to the Brody distribution two-dimensional square box in < z\ < 

function 1 and < z 2 < 1 with the potential 

XI \l (zi — Z2) 2 + S 2 similarly to the quantum 

P B (S) = (a + l)bS a exp(-bS a+1 ), 

/ + 2 \ ^ a + 1 case. The behavior of the system apparently 

^ I \a + 1/ J depends on the interaction strength param- 

which interpolates the Poisson and Wigner eter A - However, if we introduce a rescaled 

distributions. It coincides with the Pois- time T defined by r = y/\\\t, the equations 



son distribution for a = and recovers the °^ ^ e motion are reduced to 
Wigner distribution for a = 1. We use the 1 d 2 Zi _ _^ ( z i ~ Zc i) 



2 A T 2 \(za — Zo) 2 + S 2 \ 3 / 2 ' 

Brody parameter a as a measure for degree " u 1 ^ * 

1 rf 2 ^2 _ , { z 2 - z l) 

of chaotic irregularity of the system. --^ - _ ^2 + 52^3/2 ' ^ 15 J 



B. Classical dynamics 



The rescaled energy £ = i?/|A| is given by 
1 /^i\ 2 l/^2\ 2 



6 4VrfrJ + 4VdW 



Now we turn to the dynamics of the clas- 1 

± =. (16) 

sical counterpart of the two-particle system v( z i ~ z 2) 2 + $ 2 



described by the Hamiltonian (|8j) . Similarly The p i us and minus signs in Eqs.gSJ and g6|) 

to the quantum case, lengths, angular mo- correspond to positive and negative A, respec- 

mentums and masses are scaled by L, ft and tive i y> Consequently, the classical behavior of 

2m, respectively. The equations of the mo- t h e system is independent of the value of A 

tion are then given as itself. This is because in the classical system 

1 d 2 Zi X(zi — z 2 ) there is no such characteristic length as the 

2 dt {(zi — Z2) + S } / Bohr radius due to the finite Plank constant 
1 d 2 z 2 X(z 2 - zi) 



2 dt 2 {{z l -z 2 ) 2 + 8 2 Y/ 2 ' 

The A defined by A = ±2e 2 mL/h 2 is dimen- 
sionless. The total energy is represented as 

1 (dz\\ 2 1 (dz2\ 2 



(13) in the quantum system. Even if we enlarge 
the system size, we can find the equivalent 
trajectory by increasing the total energy. On 
the other hand, the quantum behavior of the 



E = + ^("J^) system depends on the value of A directly. 

A , v Now the following question arises: Even in 

+ 2 2 " ( ^ 

\J{z\ — Zc i) 2 + & 2 the quantum case, whether does the system 

The dynamics of the present system is equiv- with the same value of E but with differ- 

alent to that of a particle confined in a ent A exhibit a similar behavior to classical 
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system, especially, concerning the degree of 
chaotic irregularity of the system? Otherwise 
there is no closer correspondence between the 
classical chaos and the quantum chaos. We 
investigate this point in the present study. 

We use two kinds of the Poincare map to 
see behavior of the classical system. The first 
kind of Poincare map called Poincare map 1 
is defined in the section v 2 versus z 2 for the 
second particle taken at times when the first 
particle bounces off the left boundary of the 
well [z\ = 0). The second kind of Poincare 
map called Poincare map 2 is defined as fol- 
lows. We take coordinate / along the two 
sides of the square and a diagonal line con- 
necting (^1,^2) = (0,0) and (1,1) where the 
ridge of the potential lies (see Figj2]). / is 
normalized so that the range is < / < 1. 
The trajectory can be recorded by two val- 
ues. One is / at point where the particle is 
reflected on the hardwalls or intersects the 
line Z\ — Z2. The other is p = cos where 9 is 
the angle between the velocity vector after re- 
flection and the normal to the solid line. The 
Poincare map which records / and 6 of orbits 
reflects the property of the classical system. 

We also use MLE as a measure for degree 
of chaotic irregularity of the classical system. 
The ordinary Lyapunov exponent is defined 
as follows: we consider an orbit (p(r),q(r)) 
(denoted as the reference orbit) and a slightly 
displaced orbit from the reference orbit in the 
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FIG. 2. The coordinates of the billiard system. 

phase space. The starting point of the dis- 
placed orbit is spaced apart by a small vec- 
tor (Ap(0), Aq(0)) from (p(0),q(0)) at ini- 
tial time r = 0. The distance between the 
reference and displaced orbits is 

do = |(Ap(0),Aq(0))|. (17) 

We follow these orbits for a time interval At. 
The distance between the two orbits at r = 
At is represented as 

d 1 = |(Ap(Ar),Aq(AT))|. (18) 

Then we choose a new starting point of dis- 
placed trajectory at time r = At as 

(p(Ar),q(Ar)) 

+ ^(Ap(Ar),Aq(Ar)) (19) 

so that the distance between the new starting 
points equals d . The trajectory is followed 
up to time r = 2Ar. The new deviation of 
the displaced orbit from the reference orbit 

d 2 = |(Ap(2Ar),Aq(2Ar))| (20) 



is computed, and a second rescaled trajectory III. NUMERICAL RESULTS 
is started. This process is continued, yield- 

ing a sequence of distances d , d u d 2 , • • • . By A ' Repulsive interaction 

using these values, MLE is defined as 

We obtain energy eigenvalues by diagonal- 

1 n d 

T L = lim y^ln — , (21) izing the effective Hamiltonian rt8b, and eval- 

n^oo nAr ^ d ' 

2-1 uate the smoothed counting function N av (E) 

where n is the number of the time segment. T . T , , . T £ , . 

which gives the cumulative number ot states 

This quantity T L is, however, not suitable as , T 

below an energy b. In order to analyze 

a measure for degree of chaotic irregularity in . . £ n n n 

the structure of the level-nuctuation prop- 

the present system. Y L increases with s only 



erties, we unfold the spectrum by applying 

even because the motion of particles becomes xl nn , , r , . , . 

the well-known transformation x n = i\ av (h n ) 

faster with the increase of s, while the classi- ■> . . , . , 

to obtain a constant mean spacing, where 

n denotes the number of the energy level. 

gime as shown by the numerical results in ^ , T r i i i i ±. • xi 

i^rom the uniolded spectrum we obtain the 



cal dynamics becomes regular in high energy 
re, 

the next section. 

We introduce a novel kind of MLE defined 

as 



histogram of the NNLS distribution P(S), 
where S n = x n+ i — x n . The histogram is fit- 
ted to the Brody distribution function Pb(S) 
p/ = lim -— In — (22) in Eq.(fl2|). The integral of the Brody distri- 

2-1 bution function, 

The definition of T' T is almost the same as 

f S 

T L except that At in Eq.d2Tl) is replaced by Ib(S) = / P B {S')dS' (23) 

1 — r Jo 

a small distance ^. The definition of Y f L in 

satisfies 

Eq.([22[ is similar to Eq.Q. For Y' L we fol- 
low the two adjacent orbits while the refer- lnln[l/(l — Ib{S))\ 

ence orbit travels small distance £, and then N . ^ , , /ri . N 

s ' = (1 + a) In 5 + In 6, (24) 

evaluate d{. T' L represents the degree of expo- 



nential divergence of adjacent orbits similarly where b is given in Eq.(12). By using the 
to T L . However, Y f L depends only on geome- above relation and the least-squares fitting 
try of orbits but not on quickness of develop- method we evaluate the Brody parameter a 



ment of orbits. We employ Y f L in Eq.(22) as for the dstribution P(S) of the NNLS. Here- 
a measure of degree of chaotic irregularity in after we take S = 0.01. As an example, a 
classical mechanics. result of fitting for A = 200 is shown in Figj3| 



CO 




In S 

FIG. 3. Relation between InS and lnln[l/(l — 
7(5))] for A = 200 and 6 = 0.01. I(S) is the cu- 
mulative number of states in the unfolded spec- 
trum. The NNLS is calculated by using the 
1000th eigenstate and following 1000 eigenstates. 



The dotted line is obtained by using Eq.(24) and 
the least-squares fitting method . 




(a) level: 1000-2000, 
a = 0.221 



(d) level: 4000-5000, 
a = 0.087 
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(b) level: 2000-3000, 
a = 0.1 74 



° (e) level: 5000-6000, 4 
a = 0.014 




(c) level: 3000-4000, 
a = 0.114 



4 
(f) level: 6000-7000, 
a = 0.005 



The total energy region is divided into sev- 
eral regions. In Figj4] we show the obtained 
NNLS distribution in each region for A = 200. 
About 1000 eigenvalues are used in each re- 
gion to compute each histogram. The range 
of the used energy levels and the Brody pa- 
rameter a are shown below each panel. We 
see that the Brody parameter decreases with 
increase of the average of energy eigenvalues 
which are used to obtain the histogram. Es- 
pecially for the histograms in the panels (e) 
and (f) in Figj4] with the Brody parameter 
less than 0.015, the histograms are well fit- 
ted also by the Poisson distribution. The his- 
tograms for A = 500 is shown in Figj5| 
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FIG. 4. Histograms of NNLS for A = 200. 
Solid lines are the best fitted Brody distribu- 
tions. Each value of the Brody parameter a is 
shown below each panel. 

We see that the Brody parameter a de- 
creases with increase of the average of energy 
eigenvalues similarly to the case of A = 200, 
and moreover that the values of a are greater 
than those for A = 200 in each energy level 
region. 

We show the e— dependence of the Brody 
parameter for A = 200,500,750,1000 in 
FigJHJ For each data point we use energy 
eigenvalues in energy interval AE = 50000 




4 

(a) level: 1000-2000, 
a = 0.277 



° (d) level: 4000-5000, 4 
a = 0.2 



o 




(b) level: 2000-3000, 
a = 0.273 



(e) level: 5000-6000, 
a = 0.183 




4 4 

(c) level: 3000-4000, (f) level: 6000-7000, 

a = 0.164 a = 0.1 24 



FIG. 5. The same as Figji] except for A = 500. 

which includes about 1000 energy eigenval- 
ues. The horizontal axis denotes the av- 
eraged value of s of the used eigenstates, 

< s >. We see that the Brody parameter a 
decreases almost monotonously with increase 
of < £ >. Moreover the £— dependences of a 
is quite similar for different A especially in 

< £ >< 600. As mentioned in the previous 
section, the classical system has scaling prop- 
erty characterized by parameter e = E/\X\. 
The above results indicate that the distribu- 
tion of NNLS in quantum mechanics has the 
same scaling property on e. 



PQ 0.2 
0.1 



n 1 1 r ~ 

• X = 200 

▼ ^ = 500 

A A, = 750 
X= 1000 



V • . . 



200 400 600 800 1000 

<£> 



FIG. 6. £— dependence of the Brody parameter 
for the repulsive interaction. For each data point 
eigenvalues in the energy interval AE = 50000 
were used (Their number is about 1000). Hori- 
zontal axis denotes the average of e. 

Now we see the behavior of the classi- 
cal system with the equations of motion, 



Eq.(15). In FigsJ7j|9l we show the Poincare 



maps 1 defined in previous section for e = 
50,200, 1000, respectively, with S = 0.01. In 



Figs flQ|fl2] the Poincare maps 2 are shown. 
We have taken about 20 different initial 
points in phase space for each map. These 
Poincare maps show that the present classi- 
cal system exhibits mixed dynamics with co- 
existing KAM tri and chaotic regions. This 
is consistent with the fact that NNLS distri- 
butions in the quantum system are interme- 
diate between the Poisson and Wigner distri- 
bution. 

In Figfl3| typical trajectories ^(t)) 
are shown for (a) e = 50, (b) e = 200 and (c) 



10 




FIG. 7. Poincare map 1 for repulsive interaction FIG. 9. Poincare map 1 for repulsive interaction 
with e 50. with e 1000. 



Vi 



-2|_Zf 
-0.5 



Z 2 



0.5 




FIG. 8. Poincare map 1 for repulsive interaction FIG. 10. Poincare map 2 for repulsive interac- 

with e 200. tion with e 50. 

e = 1000, respectively. The Poincare maps enough. 

for those orbits show that they are chaotic. In order to confirm this point quantita- 



It is seen that the potential bends the trajec- tively we evaluate Y' L in Eq.(22) with £ 



tories especially near z 2 = z\ line for e = 50 0.01, which reflects irregularity of trajecto- 

and e = 200. It causes irregularity on the ries. Numerical calculation of each trajectory 

orbits. Contrastively the effect of potential is performed for time more than 1.5 xlO 5 . We 

is much less for e = 1000. The trajectory is take the average of Y f L over about 20 orbits 

composed of nearly straight lines. The orbits corresponding to the largest irregular region 

become more regular for larger e if e is large in the Poincare maps for each e. We show the 
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FIG. 11. Poincare map 2 for repulsive interac- 
tion with e — 200. 





FIG. 13. Typical trajectories {z\ (£), Z2(t)) with 
FIG. 12. Poincare map 2 for repulsive interac- ( a ) e = 50,(b) e = 200, (c) e = 1000, respec- 
tion with e = 1000. tivel Y- 



s— dependence of Y f L in Fig. 14 The decrease 
of Y f L indicates the fact that the dynamics 
becomes more regular with the increase of £, 
which is consistent with the above intuitive 
view from FigJT3| It is also consistent with 
the e— dependence of a in the quantum sys- 
tem. 

Now we calculate the ordinary MLE in 



tions for trajectories are performed for time 
more than 1.5 x 10 5 . We take the average 
of over about 20 orbits corresponding to 
the largest irregular region in the Poincare 
maps as well as for Y f L . The £— dependence 



Eq.(21) with Ar = 0.01. Numerical calcula- 



of Y L is shown in Fig{T5| We see that Y L 
increases with e. This is because motions 
of the particles become faster with the in- 
crease of £ and does not necessarily imply 
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0.45 



rv 



0.1 



2000 



maps: the area of the largest irregular region 
and the area of total region reachable for a 
particle with s. For the calculations we take 
meshes on Poincare map. A total number of 
meshes is 300 x 300. Then we count the num- 
ber of meshes which an irregular trajectory 
visits and compare it to the total number of 
meshes energetically allowed. Numerical cal- 
culations of the trajectories are performed for 



FIG. 14. e— dependence of Tr in Eq.(22). The onwin5 j j 

r l m v| — \ time more than 2.0 x 10 . Ihe e— dependence 



broken line is a guide to the eyes. 

the increase of chaotic irregularity. Therefore 
e— dependence of Y L does not directly reflect 
the degree of chaotic irregularity. 

14 




of R a for Poincare map 1 and 2 are shown in 



Figs (16] and [T7[ respectively. It is seen that 
R a increases with e for e > 100, in which 
the trajectory in two-dimensional square can 
cross the potential hill (the diagonal line). 
On the other hand, the orbits become more 



2000 



regular as seen obviously in Fig (13] when s 
increases for e > 100. Therefore R a is not a 
proper measure of irregularity in the present 
system in contrast to the other systems, in 
which R a can be adopted as a measure of ir- 
regularity P3HI5]. 



FIG. 15. The same as Fig. 14 except for in 



Eq.(21). 



We also point out that an area of the 
largest irregular region in Poincare map, 
which is adopted by several authors as a mea- 
sure of chaotic irregularity [T3l - [T5] . is irrele- 
vant for the present system. We calculate 
the ratio R a between two areas in Poincare 
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0.96 



Ra 



9 a 

it 



0.84 



500 1000 1500 2000 

£ 



FIG. 16. e— dependence of the ratio R a between 
two areas in Poincare map 1: the area of the 
largest irregular region and area of total region 
reachable for a particle with e. The broken line 
is a guide to the eyes. 



0.92 



Rc 



0.78 




FIG. 17. The same as Fig|16] except for 
Poincare map 2. 
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B. Attractive interaction 

Now we turn to results for the case of the 
attractive interaction with A < and S = 
0.01. The s— dependence of the Brody pa- 



interaction. The Poincare maps 1 are shown 



rameter is shown in Fig{T8} < s > is the aver- 
age of e of the used eigenstates. The decrease 
of the Brody parameter is seen with the in- 
crease of the scaled energy e. We see that the 
Brody parameter depends almost only on e 
and not on E and A separately, which is a sit- 
uation similar to the case of the repulsive in- 
teraction. Next, we consider the correspond- 




for £ = 5, 200, 1000 in Figs{19] - [2TJ respec- 
tively. The Poincare maps 2 are also shown 



in Figs{22]-[24} We have taken about 20 dif- 
ferent initial points in phase space for each 
map. These Poincare maps show that the 
present classical system exhibits mixed dy- 
namics with coexisting KAM tri and chaotic 
regions, as well as in the case of repulsive in- 
teraction. 




FIG. 19. Poincare map 1 for attractive interac- 
tion with e = 5. 



FIG. 18. e— dependence of the Brody pa- 
rameter for the attractive interaction. For each 
data point eigenvalues in the energy interval 
A = 50000 were used ( The number of eigen- 
values used to calculate each data point is about 
1000). Horizontal axis denotes the average of e 
of the used eigenstates. 

ing classical dynamics described by the equa- 



We calculate Y' L in Eq.Q with f = 0.01, 
which reflects the degree of chaotic irregular- 
ity. Numerical calculations of the trajectories 
are performed for time more than 1.5 x 10 5 . 



The results are shown in Figj25| where we see 
that T f L decreases with increase of e for e > 0. 
This is consistent with the fact that the dy- 
namics of the particle becomes regular when 
e increases for e > 0. On the other hand, T L 



tions of motion in Eq.(15) for the attractive in Eq.(21) increases with e as seen in Fig. 26 
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FIG. 20. Poincare map 1 for attractive interac- 
tion with e — 200. 
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-0.5 



z 2 



0.5 



FIG. 22. Poincare map 2 for attractive interac- 
tion with e = 5. 



P 



FIG. 21. Poincare map 1 for attractive interac- 
tion with e = 1000. 

where we take At = 0.01. The increase of T L 
is due to the fact that the dynamics of the 
particle becomes faster with the increase of 
s. T L does not reflect degree of chaotic irreg- 
ularity of the classical system similarly to the 
case of the repulsive interaction. 

We also calculate the ratio R a between 
two areas in Poincare maps, the area of the 
largest irregular region and the area of total 
region reachable, in the same manner used 



FIG. 23. Poincare map 2 for attractive interac- 
tion with e — 200. 

for repulsive interaction. The e— dependence 
of R a for Poincare map 1 and 2 are shown 



in Figsj27| and [28j respectively. R a increases 
with e for s > 500, while the orbit becomes 
more regular when e increases as well as in 
the case of the repulsive interaction. There- 
fore R a is not a proper measure of irregularity 
also for the case of the attractive interaction 
in present system. 
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FIG. 24. Poincare map 2 for attractive interac- 
tion with e = 1000. 



FIG. 26. The same as Fig, 25 except for MLE 



T L in Eq.(21). 




0.94 



-500 



2000 




2000 



FIG. 25. ^-dependence of T' L in Eq.(22) for the FIG. 27. ^-dependence of the ratio R a between 
attractive interaction. The broken line is a guide two areas in Poincare map 1. The broken line is 
to the eyes. a guide to the eyes. 
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FIG. 28. The same as Figj27| except for 
Poincare map 2. 
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IV. SUMMARY AND CONCLUSION 

We studied dynamics of two spinless par- 
ticles confined in a quantum wire with repul- 
sive or attractive Coulomb interaction. The 
system is reduced to a quasi-one-dimensional 
system with effective potential under the as- 
sumption that the transverse confinement is 
much stronger than the longitudinal one. 

The Coulomb interaction induces irregu- 
lar dynamics in classical mechanics. Examin- 
ing Poincare maps for the present system, we 
have found that the classical system exhibits 
mixed dynamics with coexisting KAM tori 
and chaotic regions. To see the signatures of 
quantum chaos in the corresponding quan- 
tum system we analyzed the distributions of 
the nearest neighbor level spacing (NNLS), 
which is fitted to the Brody distribution func- 
tion characterized by the Brody parameter a. 
The results indicate that they are intermedi- 
ate between the Poisson and Wigner distri- 
butions, which is consistent with the mixed 
character of the classical dynamics. 

The present classical system has a scal- 
ing property: Its dynamics is characterized 
by the rescaled energy parameter e = E/\\\, 
where A is the interaction strength parame- 
ter. Contrastingly, the quantum system has 
no such scaling property. However it has 
turned out that the distribution of NNLS in 
the quantum system has a scaling property 
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similarly to the case of classical mechanics. 
The Brody parameter a depends almost only 
on the average value of e and is insensitive to 
the value of A itself. 

In the classical system, we found that or- 
bits are more regular for larger values of s. 
The ordinary MLE T L is not suitable measure 
of chaotic irregularity for the present system, 
because they increase with e whereas the ir- 
regularity of the system decreases. We intro- 
duced a new MLE T' L , which represents a rate 
of the exponential divergence of two adjacent 
orbits (reference and displaced orbits) with 
respect to the length of the reference orbit, 
while the ordinary Lyapunov exponent de- 
scribes the one with respect to time. The de- 
pendence of T' L on e quantitatively shows the 
decrease of chaotic irregularity with increase 
of s. Therefore, the Y' L is a suitable measure 
of chaotic irregularity of the present classical 
system rather than T^. On the other hand 
in quantum system, the Brody parameter a 
decreases almost monotonously with increase 
of the average value of e in both cases of the 
repulsive and attractive interactions which 
indicates the distribution function of NNLS 
approaches to Poisson distribution with in- 
crease of e. Consequently, we have shown 
closer correspondence between the classical 
chaos and quantum chaos in the present sys- 
tem. 

We also showed that the area of the ir- 



regular region in Poincare maps are not suit- Particle-Matter Hierarchy" for its financial 

able measure of chaotic irregularity for the support. This work is partly supported by 

present system in contrast with other systems JSPS KAKENHI(Grant No. 23540459). 
in which the area has been adopted as a mea- 
sure of irregularity by several authors fT3l - fT5] . 
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